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Abstract 

We present measurements of B — > D*tv and B — > Dtv decays using 604.5 fb _1 of data 
collected at the T(4S*) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. Events are tagged by fully reconstructing one of the B mesons in hadronic modes. We 
obtain B (B+ - Wr+u) = (1.51 t° H t^9 ± 0.15)%, B (B+ - D*°t+») = (3.04 +° H ±° % ± 
0.22)%, B{B° D-t+l) = (1.01 ±g;|? ±g;if ±0.10)%, 5(5° L>*-r+^) = (2.56 ig;™ +g-3i ± 
0.10)% where the first error is statistical, the second is systematic, and the third is due to the 
uncertainty in the branching fraction for the normalization mode. 
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1. INTRODUCTION 



The semileptonic B decay to r channel, B — > D^tu, is a sensitive probe of extensions 
to the Standard Model (SM). In the SM, it occurs via an external W emission diagram with 
predicted branching fractions of (0.69 ± 0.04)% and (1.41 ± 0.07)% for the B° -> D-r + u T 
and -B — > D*~r + v T modes, respectively [l|. On the other hand, in extensions of the 
SM, such as the Two Higgs Doublet Models (2HDM) and the Minimal Supersymmetric 
Standard Model (MSSM), a charged Higgs boson (H^) may contribute to the decay 
amplitude at tree level, and the branching fraction may be modified significantly 0, y, 
0, 0, 0]- Both B — > D^tv and the purely leptonic decay B + — > t + v t have similar 
sensitivity to H^ 1 bosons with different theoretical systematics; the former suffers from 
uncertainty in the form factor, while the latter requires knowledge of the B decay constant 
/#. Therefore, the two decays provide complementary approaches to searching for H^ 1 
signatures in B decays. 

Experimentally, measurements of the B — > D^rv decays are challenging because at 
least two neutrinos are present in the final state. The Belle collaboration has previ- 
ously reported the first observation of the decay B° — > D*~t + v, by inclusively recon- 
structing the accompanying B via a 4-vector sum of all the charged and neutral tracks 
other than the D* and r daughter track candidates. The reported branching fraction is 
B(B° -»■ D*-r + u) = (2.02+°;^(stat) ± 0.37(syst))% Q. In this paper, we present new 
measurements of B° — ► D^~t + v and B + — > D^t + u decays. Here we fully reconstruct 
one of the B mesons in the event, referred to hereafter as the tag side (-Stag); and com- 
pare properties of the remaining particles, referred to as the signal side (B sig ), to those 
expected for signal and background. The method allows us to strongly suppress the com- 
binatorial backgrounds, and correctly calculate the missing mass which discriminates the 
signal from B — > D^iu background. Using a similar technique, the BaBar collaboration 
has reported the branching fractions, B(B° -> D-t + v) = (0.86 ± 0.24 ± 0.11 ± 0.06)% 
and B(B° -> D*-r + u) = (1.62 ±0.31 ±0.10 ±0.05)%, where the third uncertainty is from 
the normalization mode. They also measured distributions of the lepton momentum and 
the squared momentum transfer [si]. 

In order to avoid experimental bias, the signal region in data is blinded until the 
event selection criteria are finalized. Inclusion of the charge conjugate decays is implied 
throughout this paper. 

2. EXPERIMENT AND DATA SET 

The analysis is based on the data recorded with the Belle detector at the KEKB e + e~ 
asymmetric-energy collider operating at the center-of-mass (cm.) energy of the T(45) 
resonance. KEKB consists of a low energy ring (LER) of 3.5 GeV positrons and a high 
energy ring (HER) of 8 GeV electrons |9J]. The T(45") data set used in this analysis 
corresponds to an integrated luminosity of 605 fb -1 and contains 657 x 10 6 BB events. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a sil- 
icon vertex detector (SVD), a 50- layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintil- 
lation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a super-conducting solenoid coil that provides a 1.5 T magnetic 
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field. Muons and mesons are identified by arrays of resistive plate counters inter- 
spersed in the iron yoke (KLM). The detector is described in detail elsewhere IkJ . 

A detailed Monte Carlo (MC) simulation, based on the GEANT package [ll|, is used 
to estimate the signal detection efficiency and to study the background. Large samples of 
the signal decays are generated with the EVTGEN package 12] using the ISGW II form 
factor model [13(. To model the BB and qq(q = u, d, s, c) backgrounds, large generic BB 
and qq MC samples, corresponding to about twice the integrated luminosity of data are 
used. To further increase the BB MC statistics, we also use special BB MC samples, 
corresponding to 1.5 x 10 10 B°B° and B + B~ pairs, where events are filtered based on event 
generator information before running the time consuming GEANT detector simulation. 



3. EVENT RECONSTRUCTION AND SELECTION 
3.1. Tag-side Reconstruction 

Charged particle tracks are reconstructed from hits in the SVD and CDC. They are 
required to satisfy track quality cuts based on their impact parameters relative to the 
measured profile of the interaction point (IP) of the two beams. Charged kaons are iden- 
tified by combining information on ionization loss (dE/dx) in the CDC, Cherenkov light 
yields in the ACC and time-of-flight measured by the TOF system. For the nominal 
requirement, the kaon identification efficiency is approximately 88% and the probability 
of misidentifying a pion as a kaon is about 8%. Hadron tracks that are not identified as 
kaons are treated as pions. Tracks satisfying the lepton identification criteria, as described 
below, are removed from consideration. Electron identification is based on a combination 
of dE/dx in CDC, the response of the ACC, shower shape in the ECL, position matching 
between ECL clusters and the track, and the ratio of the energy deposited in the ECL 
to the momentum measured by the tracking system. Muon candidates are selected us- 
ing range of tracks measured in KLM and the deviation of hits from the extrapolated 
track trajectories. The lepton identification efficiencies are estimated to be about 90% 
for both electrons and muons in the momentum region above 1.2GeV/c. The hadron 
misidentification rate is measured using reconstructed Kg — > 7r + 7r~ decays and found to 
be less than 0.2% for electrons and 1.5% for muons in the same momentum region. Kg 
mesons are reconstructed using pairs of charged tracks that have an invariant mass within 
±30 MeV/c 2 of the known K s mass and a well reconstructed vertex that is displaced from 
the IP. Candidate 7's are required to have a minimum energy deposit E 1 > 50 MeV. Can- 
didate 7T° mesons are reconstructed using 77 pairs with invariant masses between 117 and 
150MeV/c 2 . For slow 7r°'s used in D* reconstruction, the minimum 7 energy requirement 
is lowered to 30 MeV. 

i?tag candidates are reconstructed in the following decay modes: B + — ► D^°ir + , 
TJ(*)° p +, JJWo a + 7JW0 jD W+ ; and B o jD (*)- 7r + ; J D(*)"p+, £>(*)- a + £>(*>-£><* )+ . Candi- 
date p + and p° mesons are reconstructed in the 7r + 7r° and tt + tt~ decay modes, by requiring 
their invariant masses to be within ±225 MeV/c 2 of the nominal p mass. We then select 
candidates by combining a p° candidate and a pion with invariant mass between 0.7 and 
1.6GeV/c 2 and require that the charged tracks form a good vertex. D meson candidates 
are reconstructed in the following decay modes: D° — > K + tt~, K + tt^tt°, K + 7t~7t + tt~ , 
K° s tt°, K° s 7i-n + , K s 7i-n + n°, K~K + , and D~ -> K+n-ji-, K + n-n-n , K° s n- , K s n-n°, 
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Kgir + 'K^'K^ , K + K~tc~ , — > K + K^it + . The D candidates are required to have 

an invariant mass Md within 4 — 5<r (a is the mass resolution) of the nominal D mass 
value depending on the mode. D* mesons are reconstructed as D* + — > D°tt + , D + tt°, 
D*° — > D°7i°, D°j, and _D* + — > -D+7. .D* candidates from modes that include a pion are 
required to have a mass difference AM = Md^ — Md within ±5MeV/c 2 of its nominal 
value. For decays with a photon, we require that the mass difference AM = M^ 7 — Mo 
be within ±20MeV/c 2 of the nominal value. 

The selection of B tag candidates is based on the beam-constrained mass Mb c = 

V -^bcam — Pb an d the energy difference AE = Eb — E\, ca , m . Here, Eb and pb are the 
reconstructed energy and momentum of the B tag candidate in the e + e~ cm. system, and 
-^beam is the beam energy in the cm. frame. The background from jet-like continuum 
events (e + e~ — > qq,q = u, d, s, c) is suppressed on the basis of event topology: we re- 
quire the normalized second Fox- Wolfram moment (R2) [HI to be smaller than 0.5, and 
I cos#th| < 0.8, where 8 t h is the angle between the thrust axis of the B candidate and that 
of the remaining tracks in the event. The latter requirement is not applied to B + — > D°tt + , 
D*°(— > D 7t )tt + and B° — > D*~(-^ D 7t~)tt + decays, where the continuum background 
is small. For the B tag candidate, we require 5.27 GeV/c 2 < Mb c < 5.29 GeV/c 2 and 
—80 MeV < AE < 60 MeV. If an event has multiple B tag candidates, we choose the 
one having the smallest x 2 based on deviations from the nominal values of AE, the D 
candidate mass, and the D* — D mass difference if applicable. The number of B + and 
B° candidates in the selected region are 1.75 x 10 6 and 1.18 x 10 6 , respectively. By 
fitting the distribution to the sum of an empirical parameterization of the background 
shape [l5| and a signal shape [l6|, we estimate that in the selected region there are 
(10.11 ± 0.03) x 10 5 (with purity=0.58) B + and (6.05 ± 0.03) x 10 5 (with purity=0.51) 
B° events, respectively. 

3.2. Signal-side Reconstruction 

In the events where a B tag is reconstructed, we search for decays of B slg into a D^*\ 
t and a neutrino. In the present analysis, the r lepton is identified in the leptonic decay 
modes, pTVv and e~Vv. We require that the charge/flavor of the r daughter particles 
and the D meson are consistent with the B slg flavor, opposite to the B tag flavor. The loss 
of signal due to B° — B° mixing is estimated by the MC simulation. 

The procedures to reconstruct charged particles (e^, /r* 1 , n"* 1 , K^) and neutral particles 
(7T°, Kg) for the signal side are the same as those used for the tag side. For 7 candidates, we 
require a minimum energy threshold of 50 MeV for the barrel, and 100 (150) MeV for the 
forward (backward) end-cap ECL. A higher threshold is used for the endcap ECL, where 
the effect of beam background is more severe. We also require that the lepton momentum 
in the laboratory frame exceeds 0.6GeV/c to ensure good lepton identification efficiency 

The decay modes used for D reconstruction are slightly different from those used for 
the tagging side: D° — > K + tt~, K + 7t~7t°, K + tt~7t + tt~ , K + 7t~7t + 7t~tt , i^^vr , K s tc~tt + , 
i^57r~7T + 7r°, and D~ — K + 7t~7t~, K + 7t~7t~tt , K s tt~ . The D candidates are required 
to have an invariant mass Md within 5cr of the nominal D mass value. D* mesons 
are reconstructed using the same decay modes as on the tagging side: D* + — > D°7r + , 
D + 7t°, and D*° — > D tc°, D 0, ~f. D* candidates are required to have a mass difference 
AM = M£) 7r ( 7 ) — Mo within 5a of the nominal value. 
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For signal selection, we use the following variables that characterize the signal decay: 
the missing mass squared in the event (M miss ), the momentum (in the cm. system) of 
the r daughter leptons (P/), and the extra energy in the ECL (Ef^). The missing mass 
squared is calculated as M miss = (E Btm - E D - E T ^ X ) 2 - {-PB tag - Pd(*) ~ Pr^x) 2 , 
using the energy and momenta of the Ptag; the Z)W candidate and the lepton from the 
r decay. The signal decay is characterized by large M^ iss due to the presence of more 
than two neutrinos in the final state. The lepton momenta (P/) distribute lower than 
those from primary B decays. The extra energy in the ECL (P^tra) i s the sum of the 
energies of photons that are not associated with either the P tag or the P S i g reconstruction. 
ECL clusters with energies greater than 50MeV in the barrel, and 100 (150) MeV in the 
forward (backward) end-cap ECL are used to calculate Ef^ a . For signal events, Ef£^ 
must be either zero or a small value arising from beam background hits, therefore, signal 
events peak at low E^^ a . On the other hand, background events are distributed toward 
higher P^ra due ^° ^ ne contribution from additional neutral clusters. We also require 
that the event has no extra charged tracks and no 7r° candidate other than those from 
the signal decay and those used in the P tag reconstruction. Table fl] summarizes the cuts 
to define the signal region. The cuts are optimized by maximizing the figure of merit 
(F.O.M.), defined as F.O.M. = N S /^N S + N B , where N S (N B ) are the number of signal 
(total background) events in the signal region, assuming the SM branching fractions for 
the Dtv and the D*tv modes. 



Cut variable B -» D°(D~)t+u B -» D*°(D*-)t + v 

= 
= 

< 1.2GeV/c 
> 1.6GeV 2 /c 4 

< 0.2 GeV 



TABLE I: Summary of the signal selection criteria. 



Number of extra tracks = 
Number of extra ir° =0 
P; <1.2GeV/c 
M 2 iss >2.0GeV 2 /c 4 
C L a < 0.2 GeV 



3.3. Signal Detection Efficiency and Expected Background 

Table [TT] lists the signal detection efficiencies, which are estimated from signal MC 
simulation, with the selection criteria shown in Table [B Taking account of the cross talks 
between B — > Dtv and B — > D*tv modes, the signal detection efficiency (e^) is defined 
as, 

Nij = 6y • Bj ■ iV tag , (1) 

where represents the yield of the generated j-th mode reconstructed in the i-th mode. 
Bj is the branching fraction of the j-th mode including the sub-decay (r and D") branch- 
ing fractions. A^ tag is the number of B events fully reconstructed on the tagging side. Ta- 
ble |TT] also shows, in parentheses, the efficiencies without cuts on M miss and Ef^ a . These 
are the two variables used to extract the signal yields. 
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T?pf*rm n mnrlp 


CrPTip'rat ,pn ttioHps 
Wt+v D*°t+v 


D°T+U 


2.55 ± 0.05 (4.87 ± 0.08) 0.90 ± 0.05 (1.75 ± 0.07) 
0.34 ± 0.01 (1.33 ± 0.02) 1.08 ± 0.03 (2.11 ± 0.04) 




Recon'd mode 


Generated modes 
D-t + u D*-t + u 


D~T + U 

D*-t+v 


3.21 ± 0.06 (6.86 ± 0.09) 0.23 ± 0.03 (0.55 ± 0.03) 
0.11 ± 0.01 (0.27 ± 0.01) 0.80 ± 0.02 (1.54 ± 0.03) 



TABLE II: Signal detection efficiency (%) matrix for B + (top) and B° (bottom) modes. The 
values in parenthesis are the efficiencies without cuts on M^ , isH and i^tra- 



According to the MC simulation, the expected number of signal (background) events 
in the signal region is 19(48) for B + -> D°r + u, 7(13) for B° -> D't+u, 18(25) for B + -> 
D*°t + v, and 7(6) for B° — > D*~t + v. The major background sources are semileptonic 
B decays, B — ► D£v, D*iv and D**iv (60-70% depending on the decay mode). The 
remaining background comes from hadronic B decays including a D meson in the final 
state. Background from qq processes are found to be small (less than one event). As 
shown in Table HT] the cross talk between B — ► Dtv and B — > D*tu arises, when a pion 
or a photon is missed in the reconstruction of D*, or when a random photon is combined 
with a D to form a fake D*. The cross-feed to other B° or B + tag samples is negligibly 
small. 

4. CALIBRATION USING THE B -> D% SAMPLE 

We use -B — > D^*>£v (£ = e//x) decays as control samples to calibrate the background 
MC simulation and to verify the Ef^ a simulation. We also use these decays to normal- 
ize the extracted signal yields. We select B — > decays using the same selection 
requirements as B — ► D(*'tv, but without the cut on the momentum of the r daughter 
lepton and with |M mis | < 1 GeV 2 /c 4 and Ef£^ < 1.2 GeV. The four calibration decay 
modes: B + -> D°£ + v, B + -> D*°£ + u, B° -> D^£+z/, and 5° -> L»*£+z/, peak around zero 
in the missing mass distributions, as shown in Figure HJ 

The yields of the calibration modes are extracted by fitting the distributions with 
expected shapes based on MC simulation for the signal and the background. The major 
background in each distribution arises from other semileptonic decays, where a pion or 
a photon is missed (i.e. B — > D*iv is reconstructed as B — » Div if the soft 7r° or 7 
from the D* is missed), or a random photon is used in D*° reconstruction (i.e. B — > Dlv 
misreconstructed as B — > D*lv). Here the two distributions for 5 + and 5° candidates 
are fitted simultaneously. Table II III lists the yields extracted for each calibration decay 
mode, which include the yields detected as cross talk; for example, the yield of D°£ + u is 
the sum of B + — ► D°£ + i> decays measured in the D°l + v and D*°£ + v distributions. When 
we compare the extracted yields with expected yields from the MC simulation, the ratio 
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of the measured to the expected yields (-R cor r) are found to be 0.75 - 0.84, depending on 
the mode. The ratios are used as scale factors to correct the normalization in the MC 
simulation for B — > D^*'£v semileptonic decays, which are the major backgrounds in the 
B — > D^tv detection. 



D°£+u D*°£+v D~£ + v D*~£+v 
Yield 1156 ± 44 2152 ± 76 338 ± 21 769 ± 35 

Efficiency (%) 8.97 ± 0.05 6.86 ± 0.08 11.3 ± 0.12 5.43 ± 0.04 

TABLE III: Yields and efficiencies of the calibration modes. 



Figure [2] compares the Ef^ a distribution for the control samples in data and the MC 
simulation after the correction. The agreement between the data and the MC simulation 
is satisfactory, and valid the Ef^ a simulation. We also confirm that the number of events 
found in the sideband of the (M miss , Ef^ a ) signal region is consistent within statistics for 
the data and the scaled MC simulation. Here the sideband is defined by E^^ a > 0.4 GeV, 
and M^ iss < 1.0GeV 2 /c 4 , for all four signal modes. 

5. SIGNAL EXTRACTION 

The B — ► Dtv and B — > D*tu signal yields are extracted using unbinned extended 
maximum likelihood fits to the two-dimensional (M^ iss , Ef^J distributions obtained after 
the selection of the signal decays. The fit components are two signal modes; B — > Dtv 
and B — > D*tu, and the backgrounds from B — > D£u, B — > D*£u and other processes. 
The likelihood is constructed as, 

-J2- N o N 

l = n e ^ ' ( 2 ) 

i=i 

where 

^(^missJ^extra) = N DrufDru{M^ iss ,Ef^+N D * TU f D * TU (M^ iss ,Ef^) 

+ N DeufD£u{M^ iss ,E^J + N D * eu f D * £u (M^ iss ,Ef^ a ) (3) 

+ Mother /other ( -^miss ; ^extra) 

Here Nj and /j(M miss , Ef^ a ) represent the number of events and the two-dimensional 
probability density function (PDF) as a function of M miss and E^ Ta , respectively, for 
process j. In the fit to the B° — > D*~t + v distribution, the Dtv cross-feed (/dtv) an d 
Div background (fniu) are not included, because their contribution are found to be small 
according to the MC simulation. The fit region is defined by (-2 GeV 2 /c 4 < Ml iss < 
8 GeV 2 /c 4 , GeV < Ef^ a < 1.2 GeV) for all four signal modes. 

The two-dimensional PDF's for D^tv and D^*'£u processes are obtained by taking 
the product of a one- dimensional PDF for each variable, as the correlations between 
M miss and i^tm are found to be small in the MC simulation. The one-dimensional 
PDF's for M^ iss are modeled by asymmetric Gaussian or double Gaussian distributions, 
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FIG. 1: Distribution of missing mass squared (M^j s ) for B + — >■ ~D°i + v (top- left), B + — >■ ~D*°l + v 
(top-right), B° —* D~£ + v (bottom-left), and B° — > D*~£ + v (bottom-right). Data are plotted 
as points with error bars, the results of the fit (solid line) along with the signal (dashed red 
line), other semileptonic decays (dotted blue line) and misidentified hadronic events (hatched 
histogram) components are also shown. 



whereas the PDF's for Ef^ a are histograms obtained from the MC simulation. The 
PDF for other background processes (/other) uses the two-dimensional histograms obtained 
from MC simulation, since correlations between the two variables are significant for these 
background processes, which mainly come from hadronic B decays. 

We fit the distributions for B° and B + tags separately. The cross talk between the 
two tags is found to be small according to the MC simulation. For each B° and B + tag, 
we then fit simultaneously the two distributions for the Drv and D*tv components. The 
ratio of the number of events found in the two distributions is constrained according to 
the efficiency matrix shown in Table [Til 
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FIG. 2: Comparison of distributions for the control samples in data and the MC simu- 

lation. 



The above procedure to extract the signal yields has been tested by "toy MC experi- 
ments": in each experiment, the number of events in each (M^ iss , Ef^^) bin is generated 
according to Poisson statistics, with the mean (//) fixed to the number of events found in 
the MC simulation, including the B — ► D^tv signals with the SM branching fraction. 
The distributions are then fit with the procedure described in the previous subsection. 
We generate 500 experiments, and we confirm that the means of the extracted yields are 
consistent with the input fi values. 

The signal extraction procedure has also been checked by performing a fit to sample 
distributions from generic MC events, which are the sum of the generic BB and qq pro- 
cesses, where semileptonic B to r decays, B — > Dru, D*ru and D**tv, are removed from 
the BB samples. For all four signal decays, the signal yields obtained are consistent with 
zero within the statistical uncertainty. 
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6. RESULTS AND SYSTEMATIC UNCERTAINTIES 



R(D°) 


= B{B + - 


^D\ + v)/B {B + - 


-> 75°£ + z/) 


R(D*°) 


= B(B + - 


^D*°t + v)/B(B + 


-> D*°£ + z/) 


R(D-) 


= B{B° - 


-> D-r + u)/B(B° - 


* ZT £ + z/) 


R(D*~) 


= B{B°- 


-> D*-r+v)/B{B° - 


-> D*-£ + u) 



In this paper, we present a relative measurement; we extract the yields of both the 
signal mode B — > D"t + i/ and the normalization mode S — > D^£ + u to deduce the four 
ratios, 

(4) 
(5) 
(6) 
(7) 

The yields of the normalization modes are extracted as described in Section [U For 
the signal modes, after finalizing the signal selection criteria and completing the studies 
in previous sections, we have opened the signal region, and performed the fits with the 
procedure described in Section[5l Figures [3] and H] show the fit results for B + — > D^tv and 
B° — > D^tis, respectively. There are excesses in the signal region for all four decay modes. 
Figure [5] shows the signal likelihood curves, while Table [IV] summarizes the results. The 
extracted yields (statistical significances) are 98.6^^(4.4), 99.8±f?;|(5.2), 17.2i£gj(2.8), 
and 25.01^(5.9), for B -> D°t + u, D*°t + u, D't+u and D*~t + v, respectively. The 
efficiency e, shown in Table IIV} corresponds to the sum of the signal yields measured in 
B -> Dtv and B -> D*tv selections. The ratio of B (B -> D^tu) to B (B -> I)W^) 
are calculated as, 



i?(-D 



N(D(*)£u) e{D(*h+v) B(t-*£vu) 



Note that the efficiency e(D^*Hu) is the average over electron and muon modes, while the 
yields are extracted for the sum of the two modes. 

Table [V] summarizes the systematic errors related to the ratio measurement, where 
reconstruction efficiency errors are largely cancel out. The following systematic errors are 
considered. 



• M^ iss shape: The systematic error due to uncertainties in the M^ iss shape is es- 
timated by varying the PDF parameters. The fitting procedure is repeated for 
each parameter variation, and relative changes in the extracted yields are added in 
quadrature. This method will give conservative estimates, as there are correlations 
in M^ iss resolutions between decay modes. 

• Ef^ a shape: The systematic error due to uncertainties in the shape is es- 
timated by varying the content of each PDF histogram bin by its ±ler statistical 
error. The fitting procedure is repeated for each bin variation, and relative changes 
in the extracted yields are added in quadrature. 

• D**lv branching fraction: The systematic errors due to uncertainties in the 
B — > D**l + v component is estimated by varying the branching fraction for each 
D** component by ±ler based on the Belle results in [l7j]. The relative change in 
the extracted yields is assigned as the systematic error. 
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M^ iss (GeV 2 /c 4 ) E extra (GeV) 



FIG. 3: Fit results for B + -> D°t+v (top) and B+ -> D*°t + u (bottom). The M^ iss (left) and 
-E'extaa (right) distributions are shown with the signal selection cut on the other variable listed 
in Table HI 



• D «-> D* cross-feed: In our nominal fitting procedure, the rates of the cross- 
feed between D and D* decays are fixed to the values in the MC simulation, for 
both the signal and normalization decays. The uncertainty is estimated by taking 
the relative change in the extracted yield for the normalization decays, when the 
cross-feed component is floated in the fit. 

• t — > tvv branching fraction: The systematic error due to uncertainties in the 
branching fraction of r decay modes is evaluated by changing the branching fractions 
by the uncertainties in the PDG values (l8j |. 

The total systematic error is the quadratic sum of all individual ones. 



10 




FIG. 4: Fit results for B° -» D~t + u (top) and B° -> D*-t + v (bottom). The M^ iss (left) and 
^eSta, (right) distributions are shown with the signal selection cut on the other variable listed 
in Table HI 



With the systematic errors shown in Table |V] the final results for the four ratios are, 

R(W) = 0.70 ±gS (9) 
R(D«) = 0.47 j£g (10) 

a(ir) = o.48 ±8:g ±g5g (ii) 

R(D*-) = 0.48 t° it ±852 , (12) 

where the first error is the statistical and the second error is the systematic. Including 
the systematic uncertainties for the yields convolved in the likelihood (Figure [5]), the 
significances of the excesses (in units of sigma) are found to be 3.8, 3.9, 2.6 and 4.7 for 
B —> D°r + v, D*°r + v, D~t + v and D*~r + v, respectively. 
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FIG. 5: Likelihood curves as a function of signal yields for B + — > D t + v (top- left) and 
B+ -> ~D*°t + v (top-right). 5° -> D~t + v (bottom-left) and 5° -> D*-t + v (bottom-right). 
Red (black) curves show the likelihood with (without) the systematic uncertainty 



Using the branching fractions for the B — > D^Hu normalization decays, reported in 
0: B{B + ^ Div) = (2.15 ± 0.22)%, B (£+ -> DHu) = (6.5 ± 0.5)%, B (B° -> £>&/) = 
(2.12 ± 0.20)%, and B (5° -> D*£z/) = (5.33 ± 0.20)%, we obtain the folowing branching 
fractions for B — > D^tu decays, 



-0.39 -0.19 

+0.69 +0.40 
-0.66 -0.47 



.15 [%] (13) 
± 0.22 [%] (14) 

±S4S ±8:S ± o.io [%] (i5) 
±S:S ±8S ± o.io [%] , (i6) 

where the first error is statistical, the second is systematic, and the third is due to the 



B{B + - 






= 1.51 


B(B + - 


-> V 


-v) 


= 3.04 


B(B° - 


-> D~V 


-v) 


= 1.01 


B(B° - 




-v) 


= 2.56 
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Quantity D"t + v D* u t + u 

N(D^t + v) 98.6±||;g 99.8±|f;§ 

e(SWr + i/) [%] 6.20 ± 0.08 3.86 ± 0.08 

pro/l 7n o +18.9 +11.0 o +10.6 +6.2 

#[/oJ 70.2 _ 18 _ 91 4b.8 _ 10 2 _ 7 2 

3.8 (4.4) 3.9 (5.2) 

B [%] 1.51 ±gj| g| ± 0.15 3.04 ±gj| ±g;g ± 0.22 

Quantity Z> _ t+i/ D*~t + v 



+7.2 
-6.3 



N(D^t + v) 17.21^ 25.0 
epWr+i/) [%] 6.86 ± 0.09 2.09 ± 0.04 

R[%] 47.6 tlH + -U 48.1 j&g 

S(S siat ) 2.6 (2.8) 4.7 (5.9) 

B[%] 1.01 ±g;g ±8:jf ± 0-10 2.56 ±||| 1°% ± 0.10 

TABLE IV: Summary of the results; extracted yields from the fitting, N, the efficiency to 
detect the signal in either of B — ► Dr(t)v and .B — > D*t(£)u selections, e, the deduced ratio 
of B (B —* D^tu) to B(B D(*'lv), R, significance of the signal with (without) systematic 
errors, S(S stat ), deducued branching fraction, B. 



Source D°r+^[%] D*°t + v[%] D~t + v^\ D*~t + u[%] 

M^ iss shape +9.10/-7.89 +9.86/-10.7 +6.39/-5.78 +5.80/-6.12 

£S shape +10.6/-7.58 +7.01/-9.73 +9.03/-7.27 +9.84/-4.97 

L>**£i/ +0.35/-0.41 +0.75/-0.02 +4.50/-2.56 +0.58/-0.28 

D ^ D* cross-feed +7.05/-6.86 +5.12/-5.34 +5.77/-6.01 +3.48/-3.37 

B(t -> Ivv) ±03 ±03 ±03 ±0.3 

Total +15.7/-12.9 +13.2/-15.4 +13.3/-11.4 +12.0/-8.58 

TABLE V: Summary of the systematic errors. 



branching fraction error of the normalization mode. 



7. SUMMARY 

Using 604.5 fb" 1 of data collected at the T(4S) resonance with the Belle detector at 
the KEKB collider, we have measured B to r semileptonic decays, by fully reconstruct- 
ing hadronic decays of the accompanying B meson. We have extracted signals for the 
four decay modes, B + -> D°r + u, B + -> D*°r + u, B° -> D't+u, and B° -> D*- T + u, 
and obtained the branching fractions listed in Section [6j The obtained branching frac- 
tions are consistent within errors with the earlier Belle result for B° — > D*~t + u [7], and 
BaBar results for the four signal modes 0. Our results are slightly higher than the SM 
expectation, however more luminosity is needed to clarify the deviation. 
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